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Abstract 

We discuss applications of the perturbative QCD approach in the exclusive non-leptonic two 
body B-meson decays. We briefly review its ingredients and some important theoretical issues 
on the factorization approach. PQCD results are compatible with present experimantal data for 
charmless B-meson decays. We predict the possibility of large direct CP asymmetry in B° — » tt + tt" 
(23 ±7%) and B° — > K + n~ (—17 ±5%). We also investigate the Branching ratios, CP asymmetry 
and isopsin symmetry breaking in radiative B — > [K* / p)j decays. 
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I. INTRODUCTION 



The aim of the study on weak decay in B-meson is two folds: (1) To determine precisely 

n □ 

the elements of Cabibbo-Kobayashi-Maskawa (CKM) matrix[l|, |2| and to explore the origin 
of CP-violation at a low energy scale, (2) To understand strong interaction physics related 
to the confinement of quarks and gluons within hadrons. 

The two tasks complement each other. An understanding of the connection between 



quarks and hadron properties is a necessary prerequeste for a precise determination o 



CKM 



matrix elements and CP-violating phases, so called Kobayashi-Maskawa(KM) phase 

The theoretical description of hadronic weak decays is difficult since nonperturbative 
QCD interactions is involved. This makes a difficult to interpret correctly data from asym- 
metric B-factories and to seek the origin of CP violation. In the case of B-meson decays 
into two light mesons, we can explain roughly branching ratios by using the factorization 
approximation 0,0]. Since B-meson is quite heavy, when it decays into two light mesons, the 
final-state mesons are moving so fast that it is difficult to exchange gluons between final-state 
mesons. So we can express the amplitude in terms of the product of weak decay constant 
and transition form factors by the factorization (color-transparancy) argument 0,1^. In this 
approach we neglect non-factorizable contributions and a power suppressed annihilation 
contributions. Because of this weakness, asymmetry of CP violation can not be predicted 
correctly. 

Recently two different QCD approaches beyond naive and general factorization assump- 
tion 0, 0, was proposed: (1) QCD-factorization in the heavy quark limit jll|, Q 
in which non-factorizable terms and cij are calculable in some cases. (2) A Novel PQCD 
approach [3, 1(| including the resummation effects of the transverse momentum carried 
by partons inside meson. In this review paper, we discuss some important theoretical issues 
in the PQCD factorization and numerical results for charmless B-decays at the section 3-7. 
In section 8 we present the PQCD results of the radiative B-decays B — > K*j, p^y. 

II. k T FACTORIZATION VS COLLINEAR FACTORIZATION 

Let's start to review shortly on the developement of theoretical methods of exclusive 
hadronic B-meson two-body decays, and make a comparision between different frameworks 
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of factorization. 

The theoretical basis how we can explain nonleptonic B-meson decays is origined from 
the color-transparancy argument 0, Q : Since b-quark decays into light quarks energetically 
(> 1 GeV), the produced quark- antiquark pair doesn't have enough time to evolve to the 
real size hadronic entity, but remains a small size bound state with a correspondingly small 
chromomagnetic moment which suppress in QCD interaction between final state mesons. 




As shown in Fig. 1, when we consider B-meson decays into two final state pions, the 
matrix elements of Oi can be expressed by a simple way as so called the naive factorization 
method 

< 7r(p2)7r(p3)|0i(^)|S(pi) > ~ < 7i(p 2 )\(d i q i ) v _ A \0 >< 7r(p 3 )|(?A)v-A|-B(>i) > 

U ® F B *{q 2 = Ml) (1) 

In this way, only factorizable part was considered, but not nonfactorizable part. In 1996, the 
generalied factorization approach was developped by A. Kamal|7j and H.-Y. Cheng which 
included non-factorizable contributions into the effective wilson coefficients by assuming 
NF = x® F. The generalized factorization has a weak point to predict the CP asymmetry, 
since they considered non-factorizable part to be real, however it is complex in general. After 
then, QCD-factorization was proposed by Beneke et al.jllj], which is an improved form of 



naive factorization approach. When we consider B — > M\M.2 with recoiled M\ and emitted 
M2 (light or quarkonium) , soft gluon exchanged effects are confined to [BM\) system and 
only hard interactions between [BM\] and M 2 survive in mj, — > 00 limit which is calculable 
perturbatively. The decay amplitude can be written as: 

Decay Amp = Amp ( „ fact) ® [1 + 0(a s ) + 0(^£)]. (2) 

In principle, nonfactorizable part and a« are calculable in the heavy quark limit within the 
leading twist, however it is diverged at the end point with the twist-3 contributions and 
even the leading twist contribution in the annihilation diagram. To solve this end-point 
sigularity problem, PQCD approach was proposed by Keum et al. 13j], in which hard gluon 
exchanged contributions is dominant even in the {BM\) transition form factor. We will 
discuss the detail of PQCD approach in the next section. 



Now we explain how to derive collinear and hr factorization theorems for the pion form 
factor involved in the scattering process 7r(-Pi)7*((?) — > ^(Pz)- The momenta are chosen 
in the light-cone coordinates as Pi = (P* ,0,0^), Pi = (0,P 2 ~,Ot), and Q 2 = —q 2 . At 
leading order, 0(a s ), shown in Fig. 2(a), the hard kernel is proportional to if ^(xi, 22) oc 
— 1/ {x\P\ — X2P2) 2 = V (xiX2Q 2 )- Here x\ and Xi are the parton momentum fractions carried 
by the lower quarks in the incoming and outgoing pions, respectively. At next-to-leading 
order, 0(a 2 ), collinear divergences are generated in loop integrals, and need to be factorized 
into the pion wave function. In the collinear region with the loop momentum I parallel 
to Pi, we have an on-shell gluon I 2 ~ P 2 ~ 0(A 2 ) with the hierachy of the components, 

1+ ~ p+ > i T ^ a > r ~ A 2 /pf. 




FIG. 2: (a) Lowest-order diagram for F n . (b) Radiative correction to (a). 
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An example of next-to- leading-order diagrams is shown in Fig. 2(b). The factorization of 
Fig. 2(b) is trivial: one performs the Fierz transformation to separate the fermion flows, so 
that the right-hand side of the cut corresponds to the lowest-order hard kernel H™. Since 
the loop momentum / flows into the hard gluon, we have the gluon momentum X\P\ — X2P2 + I 
and 

ff(0) (x — (3) 

(xiPi - x 2 P 2 ) 2 + 2xiP 1 + Z- - 2x 2 P 2 ~l + + 21+1- -I 2 ' v ' 

Dropping l~ and It as a collinear approximation, the above expression reduces to 

^ (0)(6 ' X2) * 2x 1 x 2 P 1 + P 2 + 2x 2 P 2 l + = I^Q 2 ' W 
where 6 = x\ + l + / P\ is the parton momentum fraction modified by the collinear gluon 
exchange. The left-hand side of the cut then contributes to the 0(a s ) distribution amplitude 
<f>Tp(£i), which contains the integration over l~ and It- Therefore, factorization to all orders 
gives a convolution only in the longitudinal components of parton momentum, 

F w = J ^1^2^(6)^(6,6)^(6) • (5) 

In the region with small parton momentum fractions, the hard scale X1X2Q 2 is not large, 
n this case one may drop only l~ , and keep It in H^°\ This weaker approximation gives 

which acquires a dependence on a transverse momentum. We factorize the left-hand side of 
the cut in Fig. 2(b) into the 0(a s ) wave function 0i (6>^t)j which involves the integration 
over l~ . It is understood that the collinear gluon exchange not only modifies the momentum 
fraction, but introduces the transverse momentum dependence of the pion wave function. 
Extending the above procedure to all orders, we derive the fcy factorization, 

F-k = J d^ 1 d^ 2 d 2 k 1T d 2 k2T(pA^i, hr)H (f 1, 6, hr, fer)07r(6> hr) ■ (7) 
III. INGREDIENTS OF k T FACTORIZATION APPROACH 



Factorization in PQCD: The idea of pertubative QCD is as follows: When heavy B- 
meson decays into two light mesons, the hard process is dominant. Since two light mesons fly 
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so fast with large momentum, it is reasonable assumptions that the final-state interaction is 
not important for charmless B-decays. Hard gluons are needed to boost the resting spectator 
quark to get large momentum and finally to hadronize a fast moving final meson. So the 
dominant process is that one hard gluon is exchanged between specator quark and other 
four quarks. 

Let's start with the lowest-order diagram of B — > Kit. The soft divergences in the 
B — > 7i form factor can be factorized into a light-cone B meson wave function, and the 
collinear divergences can be absorbed into a pion distribution amplitude. The finite pieces 
of them is absorbed into the hard part. Then in the natural way we can factorize amplitude 
into two pieces: G = H(Q,fi) ® $(m, /i) where H stands for hard part which is calculable 
with a perturbative way. $ represents a product of wave functions which contains all the 
nonperturbative dynamics. 

PQCD adopt the three scale factorization theorem based on the perturbative QCD 
formalism by Brodsky and Lepage [3], and Botts and Sterman with the inclusion of 
the transverse momentum components carried by partons inside meson. 

We have three different scales: electroweak scale: Mw, hard interaction scale: t ~ 
0(y(Arn.&)), and the factorization scale: 1/6 where b is the conjugate variable of parton 
transverse momenta. The dynamics below 1/6 is completely non-perturbative and can be 
parameterized into meson wave funtions which are universal and process independent. In 



our analysis we use the results of light-cone distribution amplitudes (LCDAs) by Ball 20L 12 1 1 
with the light-cone sum rule calculation. 
The ampltitude in PQCD is expressed as 



(M 1 M 2 \C k (t)O k \B) = J [dx] J 



<Pb 



®*(x 2 ,b 2 )$* M (x 3 ,b 3 )C k (t) 



® H k {{x], {6}, M b )*b{x u h) S t ({x}) e ~s(MM,M B ) (8) 

with the sudakov suppressed factor: S = SgixiP^, b\) + Sm x {x^P^ , b 2 ) + Smi((1 — 
X2)P-2 r , ^2) + ••• and the threshold resummation factor S t (x). Here C(t) are Wilson coef- 
ficients, are meson LCDAs and variable t is the factorized scale in hard part. 
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FIG. 3: The diagrams generate double logarithm corrections for the sudakov resummation. 



Sudakov Suppression Effects: When we include k±, the double logarithms ln 2 (P6) 
are generated from the overlap of collinear and soft divergence in radiative corrections to 
meson wave functions(See figure 3), where P is the dominant light-cone component of a 
meson momentum. The resummation of these double logarithms leads to a Sudakov form 
factor exp[— s(P, b)] in Eq.(8), which suppresses the long distance contributions in the large 
b region, and vanishes as b > 1/Aqcd- This suppression renders k\ flowing into the hard 



exp[-s(Q,i>)] 





FIG. 4: (a)Sudakov suppression factor (b)Fractional contribution to the B — > tt transition form 
factor F function of a s (t)/n. 

amplitudes of order 



< kl >r 



0(AM 



B 



(9) 



The off-shellness of internal particles then remain of O(AMg) even in the end-point region, 
and the singularities are removed. This mechanism is so-called Sudakov suppression (See 
figure 4-a). 

Du et al. have studied the Sudakov effects in the evaluation of nonfactorizable amplitudes 



22j | . If equating these amplitudes with Sudakov suppression included to the parametrization 



lat the corresponding cutoffs are located in the reasonable range 



in QCDF, it was observed t 
proposed by Beneke et al. |12|. Sachrajda et al. have expressed an opposite opinion on the 
effect of Sudakov suppression in j^J. However, their conclusion was drawn based on a very 
sharp B meson wave function, which is not favored by experimental data. 

Here I would like to comment on the negative opinions on the large k\ ~ 0(AM B ). It 
is easy to understand the increase of k\ from 0(A 2 ), carried by the valence quarks which 



just come out of the initial meson wave functions, to 0(AMb), carried by the quarks which 
are involved in the hard weak decays. Consider the simple deeply inelastic scattering of 
a hadron. The transverse momentum k± carried by a parton, which just come out of the 
hadron distribution function, is initially small. After infinite many gluon radiations, kj_ 
becomes of 0(Q), when the parton is scattered by the highly virtual photon, where Q is 
the large momentum transfer from the photon. The evolution of the hadron distribution 
function from the low scale t o Q is described by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equation 2J,|25[. The mechanism of the DGLAP evolution in DIS is similar 
to that of the Sudakov evolution in exclusive B meson decays. The difference is only that 
the former is the consequence of the single-logarithm resummation, while the latter is the 
consequence of the double-logarithm resummation. 

By including Sudakov effects, all contributions of the B — ► 7r form factor comes from the 
region with a s jn < 0.3 [15] as shown in Figure 4(b). It indicate that our PQCD results are 
well within the perturbative region. 



FIG. 5: The diagrams generate double logarithm corrections for the threshold resummation. 



Threshold Resummation: The other double logarithm is a,ln 2 (l/a;) from the end 

n 

point region of the momentum fraction x |26lj . This double logarithm is generated by the 
corrections of the hard part in Figure 5. This double logarithm can be factored out of the 
hard amplitude systematically, and its resummation introduces a Sudakov factor S t {x) = 
1.78[x(l — x)] c with c = 0.3 into PQCD factorization formula. The Sudakov factor from 
threshold resummation is universal, independent of flavors of internal quarks, the twists and 
topologies of hard amplitudes, and the decay modes. 



27 



281 ] arise from different subpro- 



Threshold resummation 1261 and k± resummation 
cesses in PQCD factorization and suppresses the end-point contributions, making PQCD 
evaluation of exclusive B meson decays reliable. Without these resummation effects, the 
PQCD predictions for the B — > K form factors are infrared divergent. If including only 
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k± resummation, the PQCD predictions are finite. However, the two-parton twist-3 contri- 
butions are still huge, so that the B — > K form factors have an unreasonably large value 
F BK ~ 0.57 at maximal recoil. The reason is that the double logarithms a s ln 2 x have not 
been organized. If including both resummations, we obtain the reasonable result F BK ~ 0.35 
as shown in Figure 6.. These studies indicate the importance of resummations in PQCD 
analyses of B meson decays. In conclusion, if the PQCD analysis of the heavy-to-light form 
factors is performed self-consistently, there exist no end-point singularities, and both twist- 2 
and twist-3 contributions are well-behaved. 



Amplitudes 


twist-2 contribution Twist-3 contribution 


Total 


Re{f,F T ) 
Im(f n F T ) 


3.44 • 10~ 2 5.00 • 10~ 2 


8.44 • 10~ 2 


Re(UF p ) 
Im(UF p ) 


-1.26 • 10" 3 -4.76 • 10~ 3 


-6.02 • 10~ 3 


Re{f B F?) 
Im{f B F p ) 


2.52 • 10~ 6 -3.30 • 10" 4 
8.72 • 10~ 7 3.81 • 10~ 3 


-3.33 • 10~ 4 
3.81 • 10~ 3 


Re(M T ) 
Im(M T ) 


7.26 • 10~ 4 -1.39 • 10~ 6 
-1.62 • 10~ 3 -2.91 • 10~ 4 


-7.25 • 10~ 4 
1.33 • 10~ 3 


Re(M p ) 
Im(M p ) 


-1.67 -lO" 5 -1.47 -lO" 7 
-3.52 • 10~ 5 6.56 • 10~ 6 


1.66 • 10~ 5 
-2.87 • 10~ 5 


Re{M p ) 
Im(M p ) 


-7.37 • 10~ 5 2.50 • 10~ 6 
-3.13 • 10~ 5 -2.04 • 10" 5 


-7.12 • 10~ 5 
-5.17 • 10~ 5 



TABLE I: Amplitudes for the — > it + it decay where F (M) denotes factorizable (nonfactoriz- 
able) contributions, P (T) denotes the penguin (tree) contributions, and a denotes the annihilation 
contributions. Here we adopted = 80°, R b = 0.38, = 1.4 GeV and u B = 0.40 GeV. 



Power Counting Rule in PQCD: The power behaviors of various topologies of dia- 
grams for two-body nonleptonic B meson decays with the Sudakov effects taken into account 
has been discussed in details in |29l The relative importance is summarized below: 



2?77q A. 

emission : annihilation : nonfactorizable = 1 : : , (10) 

M B M B ' K ' 
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J 0.1 

X 



0.6 



0.4 



0.2 



— Sudakov, F^ >K (0)=0.35 

— No Su. in b 1 and b 2 , 1.16 
...No Su. in b r 0.36 
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0.3 



T-t-T 
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FIG. 6: Sudakov suppression and threshold resummation effects in B — ► K transition form factor 



with m being the chiral symmetry breaking scale. The scale m appears because the anni- 
hilation contributions are dominated by those from the (V — A)(V + A) penguin operators, 
which survive under helicity suppression. In the heavy quark limit the annihilation and non- 
factorizable amplitudes are indeed power-suppressed compared to the factorizable emission 
ones. Therefore, the PQCD formalism for two-body charmless nonleptonic B meson decays 
coincides with the factorization approach as Mb — > oo. However, for the physical value 
Mb ~ 5 GeV, the annihilation contributions are essential. In Table 1 and 2 we can easily 
check the relative size of the different topology in Eq. (JlOJ) by the peguin contribution for W- 
emission (f n F p ), annihilation(/^i 71 ( f ) and non-factorizable(M p ) contributions as shown in 
Figure 7. Specially we show the relative size of the different twisted light-cone-distribution- 
amplitudes (LCDAs) for each topology. We have more sizable twist-3 contributions in the 
factorizable diagram. 

Note that all the above topologies are of the same order in a s in PQCD. The nonfactor- 
izable amplitudes are down by a power of because of the cancellation between a pair 

of nonfactorizable diagrams, though each of them is of the same power as the factorizable 
one. I emphasize that it is more appropriate to include the nonfactorizable contributions 
in a complete formalism. The factorizable internal-VF emisson contributions are strongly 
suppressed by the vanishing Wilson coefficient a? in the B J/ipK^ decays jsij, so that 
nonfactorizable contributions become dominant j3l| . In the B — ► Dn decays, there is no soft 
cancellation between a pair of nonfactorizable diagrams, and nonfactorizable contributions 
are significant 3cl ]. 
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FIG. 7: Feynman diagrams for B — > rnr and i^7T. 



In QCDF the factorizable and nonfactorizable amplitudes are of the same power in l/m&, 
but the latter is of next-to-leading order in a s compared to the former. Hence, QCDF 
approaches FA in the heavy quark limit in the sense of a s — > 0. Briefly speaking, QCDF 
and PQCD have different counting rules both in a s and in l/m&. The former approaches 
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FA logarithmically (a s oc 1/ln 777.5 - > 0), while the latter does linearly (l/m& — * 0). 



Amplitudes 


Left-handed gluon exchange Right-handed gluon exchange 


lotal 


Re(f K F T ) 

^ t tpT\ 
lm(fKr ) 


7.07 • 1(T 2 3.16 • 1(T 2 


1.02 • 10" 1 


Re{f K F p ) 
lm{fKr ) 


-5.52 • 1(T 3 -2.44 • 1(T 3 


-7.96 • 10~ 3 


J nm / J* T7'P\ 

Im(f B F a ) 


4.13 • 10~ 4 -6.51 • 10~ 4 

70 1 n — 3 1 ^0 1 n — 3 
i. to ■ 1U l.Do • ID 


-2.38 • 10~ 4 

A A 1 1 n — 3 
4.41 • 1U 


Re(M T ) 

/ i\/rT\ 

Imylvl ) 


7.06 • 10~ 3 -7.17 • 10~ 3 

-1.1U * 1U l.oO • 1U 


-1.11 • 10~ 4 
Ka 1 n — 3 


Re(M p ) 
Im(M p ) 


-3.05 • 10~ 4 3.07 • 10~ 4 
4.50 • 10~ 4 -5.29 • 10~ 4 


2.17- 10~ 6 
-7.92 • 10~ 5 


Re(M p ) 
Im(M p ) 


2.03 • 10~ 5 -1.37 -10" 4 
-1.45 -lO" 5 -1.27 -10" 4 


-1.16 • 10~ 4 
-1.42 • 10~ 4 



TABLE II: Amplitudes for the — > K + -k~ decay where F (M) denotes factorizable (nonfactoriz- 
able) contributions, P (T) denotes the penguin (tree) contributions, and a denotes the annihilation 
contributions. Here we adopted 03 = 80°, = 0.38. 



IV. IMPORTANT THEORETICAL ISSUES 



End Point Singularity and Form Factors: If calcula ting the B — ► 71 form factor 



F at large recoil using the Brodsky-Lepage formalism ], a difficulty immediately 

occurs. The lowest-order diagram for the hard amplitude is proportional to l/(xxx|), x\ 
being the momentum fraction associated with the spectator quark on the B meson side. If 
the pion distribution amplitude vanishes like X3 as £3 — > (in the leading-twist, i.e., twist- 2 
case), F BlT is logarithmically divergent. If the pion distribution amplitude is a constant as 
x 3 — > (in the next-to- leading- twist, i.e., twist-3 case), even becomes linearly divergent. 
These end-point singularities have also appeared in the evaluation of the nonfactorizable and 
annihilation amplitudes in QCDF mentioned above. 
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When we include small parton transverse momenta k±, we have 
1 1 



(x 3 M| + kj ± ) [ Xl x 3 M 2 B + (k 1± - k 3± ) 



xi x\M B 

and the end-point singularity is smeared out. 

In PQCD, we can calculate analytically space-like form 
and also time-like form factors for the annihilation process 



'actors for B 

3Q- 



in 



P, V transition 



Strong Phases: While stong phases in FA and QCDF come from the Bander-Silverman- 
Soni (BSS) mechanism 34] and from the final state interaction (FS D, t he dominant strong 
phase in PQCD come from the factorizable annihilation diagram 3,130] (See Figure 8). 
In fact, the two sources of strong phases in the FA and QCDF approaches are strongly sup- 
pressed by the charm mass threshold and by the end-point behavior of meson wave functions. 
So the strong phase in QCDF is almost zero without soft-annihilation contributions. 



Vertex diagram in QCDF Annihilation diagram in PQCD 

FIG. 8: Different sourses of strong phase: (a)BSS mechanism, (b) Final State Interaction, and (c) 
Factorizable annihilation. 




BSS — mechanism 



Dynamical Penguin Enhancement vs Chiral Enhancement: As explained before, 
the hard scale is about 1.5 GeV. Since the RG evolution of the Wilson coefficients C^(t) 
increase drastically as t < Mb/2, while that of C\^{t) remain almost constant as shown in 
Figure 9, we can get a large enhancement effects from both wilson coeffi cents and matrix 
elements in PQCD. 

In general the amplitude can be expressed as 

Amp ~ [ai,2 ± 04 ± m^ ([i)a§\ ■ < Kn\0\B > (12) 

with the chiral factors m^(/i) = m 2 p /[mi(fi) + m2(fi)} for pseudoscalr meson and m$ = my 
for vector mesons. To accommodate the B — + Kir data in the factorization and QCD- 
factorization approaches, one relies on the chiral enhancement by increasing the mass m Q to 
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FIG. 9: Dynamical enhancement of Wilson coefficents en (i=l,4,6). 



as large values about 3 GeV at [i = scale. So two methods accomodate large branching 
ratios of B — > Kir and it is difficult for us to distinguish two different methods in B — > PP 
decays. However we can do it in B — > PV because there is no chiral factor in LCDAs of 
the vector meson. We can test whether dynafrrffcal enhancement or chiral enhancement is 
responsible for the large B — > Kir branching ratios by measuring the B — > <f)K modes. In 
these modes penguin contributions dominate, such that their branching ratios are insensitive 
to the variation of the unitarity angle According to recent works by Cheng at al. j^, 
the branching ratio of B — > <pK is (2 — 7) x 10~ 6 including 30% annihilation contributions in 
the Q CD-factorization approach (QCDF). However PQCD predicts 10 x 10" 6 13, Q- For 
B (pK* decays, QCDF gets about 9 x 10~ 6 [36], but PQCD have 15 x lO ^j^J Because 
of these small branching ratios for B — > PV and VV decays in the QCD-factorization 
approach, they can not globally fit the experimental data for B — > PP, VP and VV modes 
simultaneously with same sets of free parameters (ph,4>h) and (pa,4>a) |^Z|- 



Fat Imaginary Penguin in Annihilation: There is a falklore that the annihilation 
contribution is negligible compared to W-emission one. In this reason the annihilation 
contribution was not included in the general factorization approach and the first paper 
on QCD-factorization by Beneke et al. [JjJ. In fact there is a suppression effect for the 
operators with structure (V — A)(V — A) because of a mechanism similar to the helicity 
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suppression for ir — > [w^. However annihilation from the operators 05,6,7,8 with the structure 
(S — P)(S + P) via Fiertz transformation survive under the helicity suppression and can get 
large imaginary value. The real part of factorized annihilation contribution becomes small 
because there is a cancellation between left-handed gluon exchanged one and right-handed 
gluon exchanged one as shown in Table 2. This mostly pure imaginary value of annihilation 
is a main source of large CP asymmetry in B — > 7r + 7r~ and K + n~ . In Table 6 we summarize 
the CP asymmetry in B — > K(7r)7r decays. 



V. NUMERICAL RESULTS 



Branching ratios in Charmless B-decays: The PQCD approach allows us to cal- 
culate the amplitudes for charmless B-meson decays in terms of ligh-cone distribution 
amplitudes upto twist-3. We focus on decays whose branching ratios have already been 
measured. We take allowed ranges of shape parameter for the B-meson wave funtion as 
ujb = 0.36 — 0.44 which accomodate to reasonable form factors, F B7T (0) = 0.27 — 0.33 and 
F BK (0) = 0.31-0.40. We use values of chiral factor with m£ = 1.3GeV and m$ = 1.7GeV. 
Finally we obtain branching ratios for B K(-k)tt 0,0, K( P E> 42 1 ^*</>Q and iCvrQ, 
which is well agreed with present experimental data in Table 3-5. 

CP Asymmetry of B — > tttt, Ktt: Because we have a large imaginary contribution from 
factorized annihilation diagrams in PQCD approach, we predict large CP asymmetry (~ 25%) in 
B° — > 7r + 7r _ decays and about —15% CP violation effects in B° — > K + ir~. The detail prediction 
is given in Table 6. The CP asymmetry is defined as followings: 

N(B -> /) - N(B -» /) 



A C p{At) 



N(B -> /") + N(B -» /) 

5/ sin(Am d At) - C f Cos(Am d At). (13) 



Here we notice that the relation between two different definitions: A f (Belle) = —Cf(BaBar). In 
our analysis we used the Belle notation. The precise measurement of direct CP asymmetry (both 
magnitude and sign) is a crucial way to test factorization models which have different sources of 
strong phases. Our predictions for CP-asymmetry on B — > K(tt)tt have a totally opposite sign 
to those of QCD factorization. Recently it was confirmed as the first evidence of the direct CP- 
violation in B-decays that the DCP asymmetry in B — ► K^tt^ decay are —10.1 ± 2.6% with 3.9a 
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TABLE III: Branching ratios of B -> vrvr, K?rand if-fiT decays with <fo = 80°, i? 6 = yV + r/ 2 = 0.38. 
Here we adopted ml = 1.3 GeV, = 1.7 GeV and 0.36 < uj b < 0.44. Unit is 10~ 6 . 



Modes 


CLEO BELLE BABAR 


World Av. 


PQCD 


7T + 7r _ 
„+„0 

7r°vr 


4 - 5 ±L2-o:4 4.4 ±0.6 ±0.3 4.7 ±0.6 ±0.2 

A rr - h 1 .8— |— . D cr O l 1 o l n c C O 1 r\ G 1 n A 

4.5_ x 6 Iq 7 5.3 ±1.3 ±0.5 5.8 ± 0.6 ± 0.4 
< 4.4 2.32toittoil 1.7 ±0.32 ±0.10 


4.55 ± 0.44 
5.20 ± 0.79 
2.01 ± 0.43 


5.93 - 10.99 
2.7 2 — 4.79 
0.1 - 0.65 


K°7T° 


18.0±| f±J;| 18.5 ± 1.0 ± 0.7 17.9 ± 0.9 ± 0.7 

1 Q Q+3. 7+2.1 r)r) r\ i i n j_ 1 1 n J- i q _i_ 1 n 

io.o_g g zz.u =t i.y =t i.i zo.u zt i.o ± i.u 
12.91^};? 12.8 ±1.4+};^ 12.0 ± 0.7 ± 0.6 
12.8lfjjl};][ 12.6 ± 2.4 ± 1.4 11.4 ± 0.9 ± 0.6 


18.2 ±0.8 

ZZ.o zt 1.4 

12.6 ± 1.1 

12.3 ± 1.7 


12.67 - 19.30 

14. 4o — ZD. ZD 

7.87 - 14.21 
4.46 - 8.06 




< 0.8 < 0.7 < 0.6 

< 3.3 <3.4 1.45 ±0.50 ±0.11 

< 3.3 <3.2 1.19 ±0.38 ±0.13 


< 0.6 
1.45 ±0.50 ±0.11 
1.19 ±0.38 ±0.13 


0.06 
1.4 
1.4 



TABLE IV: Branching ratios of B -> cj)K ^and K*vr decays with 3 = 80°, i? 6 = yV + r/ 2 = 0.38. 
Here we adopted = 1.3 GeV and = 1.7 GeV. Unit is 10 -6 . 



Modes 


CLEO 


BELLE 


BABAR 


World Av. 


PQCD 


0K ± 


5.5t 2 ;g±0.6 9.4 ±1.1 ±0.7 


10.0ig| ± 0.5 


9.3 ±0.8 


8.1 


- 14.1 


<j)K° 


5.4+| 7 ±0.7 9.0 ±2.2 ±0.7 


7.6±J;| ± 0.5 


7.7 ± 1.1 


7.6 


- 13.3 


fK** 
4>K*° 


10 6 +6 ' 4+LS 
iu -°-4.9-1.6 

11 k+4.5+1.8 
li -' J -3.7-1.7 


n 7 2.1+0.7 
-1.9-1.0 

10 n +L6+0 - 7 

iu - u -l. 5-0.8 


12.1+I 1 ± 1.1 
ll.l±i;| ± 0.8 


9.4 ±1.6 
10.7 ± 1.1 


12.6 
11.5 


- 21.2 

- 19.8 


K*°TT ± 


7.6jf;{j ± 1.6 
16tt ± 2 


10 4+4.2+4.1 
iy -^-3. 9-7.1 

< 30 


15.5 ±3.4 ±1.8 


12.3 ±2.6 
16 ±6 


10.2 
8.0 


- 14.6 

- 11.6 




< 31 






< 31 


2.0 


- 5.1 


K*°ir° 


< 3.6 


< 7 




< 3.6 


1.8 


-4.4 



deviations from zero in Belle Coll., and —13.3 ±3.1% with 4.2cr, which is in a good agreement with 
PQCD result 3- 
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VI. EXTRACTION OF 2 (= a) FROM B -> Tr+Tr" 

Even though isospin analysis of B — > 7T7r can provide a clean way to determine 02 , it might be 
difficult in practice because of the large uncertainty of the branching ratio of B° — > 7r°7r°. In reality 
in order to determine 02, we can use the time-dependent rate of B°(t) — > 7r + 7r~. Since penguin 
contributions are sizable about 20-30 % of the total amplitude, we expect that direct CP violation 
can be large if strong phases are different in the tree and penguin diagrams. 

In our analysis we use the c-convention. The ratio between penguin and tree amplitudes is 
Rc = \Pc/T c \ and the strong phase difference between penguin and tree amplitudes 5 = 5p — 5t- 
The time-dependent asymmetry measurement provides two equations for C n7T and S nn in terms 
of three unknown variables R C ,S and 02|49J]. Since PQCD provides us R c = 0.23^q'q 5 and 
—41° < 5 < —32°, the allowed range of 02 at present stage is determined as 55° < 02 < 100° as 
shown in Figure IVTI 

According to the power counting rule in the PQCD approach, the factorizable annihilation 
contribution with large imaginary part becomes subdominant and give a negative strong phase 
from —ii\:5{k' 2 1 _ — xMg). Therefore we have a relatively large strong phase in contrast to the 
QCD-factorization (8 ~ 0°) and predict large direct CP violation effect in B° — > ir + ir~ with 
A cp (B° — ► 7r + 7r~) = (23 ± 7)%, which will be tested by more precise experimental measurement 



Quatity 


Experiment 


PQCD 


QCDF[45j 


Br(n + iT~ ) 
Br(n ± K~^) 

Br(vr±KT) 
2Br(-K K°) 

2 Br(7r°K' ± ) 
Br(TT±K°) 

t(B+) Br(n^K ± ) 
r(B u ) Br(7r±i4' <J ) 


0.25 ± 0.04 
1.05 ±0.27 
1.25 ±0.22 
1.07 ±0.14 


0.30 - 0.69 
0.78 - 1.05 
0.77-1.60 
0.70 - 1.45 


0.5 - 1.9 
0.9 - 1.4 
0.9 - 1.3 
0.6 - 1.0 



TABLE V: Ratios of CP-averaged rates in B — ► Kir, nir decays with 03 = 80°, Rb = 0.38. Here we 
adopted = 1.3 GeV and = 1.7 GeV. 
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Direct A CP {%) 


BELLE 


BABAR 


PQCD 


QCDF 


7T + 7T~ 


58 ± 15 ± 7 


9± 15 ±4 


16.0 ~ 30.0 


-6 ±12 


7T+7T 


-14 ± 24+1 


1 ± 10 ± 2 


0.0 


0.0 


vr vr 


43 ± 51+^ 


12 ± 56 ± 6 


20.0 ~ 40.0 





7T+K- 


-10.1 ±2.5 ±0.5 


-13.3 ±3.0 ±0.9 


-12.9 21.9 


5 ±9 


7T°K- 


4±5±2 


6.0 ±6.0 ±1.0 


-10.0 ~ -17.3 


7±9 




7+9+1 
'-8-3 


-8.7 ± 4.6 ± 1.0 


-0.6 1.5 


1 ± 1 


ir°K° 


16 ± 29 ± 5 


-6 ±18 ±6 


-0.90 1.03 


-3.6 ~ 0.8 



TABLE VI: CP-asymmetry in B -> Kit, -kit decays with (j) 3 = 40° ~ 90°, R b = V ' P 2 + V 2 = 0.38. 
Here we adopted = 1.3 GeV and uiq =1.7 GeV. 

2 i i 



Belle: 


A=0.58+ 0.17, 


S=-1.00 ± 0.22 


BaBar: 


A=0.09+ 0.16. 


S= -0.30t 0.17 




FIG. 10: Plot of A W7r versus S n7r for various values of <f>2 with <f>\ = 24.3°, 0.18 < R c < 0.30 and 
-41° < 8 < -32° in the pQCD method. 



within two years. 

In the numerical analysis, though the data by Belle collaboration! 541 is located ourside allowed 

nn 

physical regions, we considered the averaged value of recent measurements |54. 15q : 

• S nn = -0.30 ± 0.17 ± 0.03 (BaBar), S n7T = -1.00 ± 0.21 ± 0.07 (Belle); 

• A nn = 0.09 ± 0.15 ± 0.04 (BaBar), A nn = 0.58 ± 0.15 ± 0.07 (Belle). 

The central point of averaged data corresponds to fa = 78° in the PQCD method. Even if the data 
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n 

by Belle collaboration [5J] is located ourside allowed physical regions, we can have allowed ranges 
with 2 a bounds, but large negative S and R c > 0.4 is prefered 



VII. EXTRACTION OF 3 (= 7) FROM B° -> K+vr- AND -> i^°vr + 

By using tree-penguin interference in B° — ► l^ + 7r~(~ T + P ) versus -B + — » i^°7r + (~ 
CP-averaged 5 — > X7r branching fraction may lead to non-trivial constaints on the <f>% angle 
In order to determine 03, we need one more useful information on CP- violating rate differences 

Let's introduce the following observables : 

Br(B° K+ir-) t + 2 
#K = p-, p , 77n TT = 1-2 tk cosd cosfo + r K 

> sm 2 3 (14) 

r( J B° -> - r(s° K+7T-) 

^ cp (i? — > K + 7T _ ) i?x = -2r K sinfo sind. (15) 



A. 



where rx = \T jP | is the ratio of tree to penguin amplitudes and 5 = 5t> — 5 pi is the strong phase 
difference between tree and penguin amplitides. After eliminate sin5 in Eq.(8)-(9), we have 

Rk = 1 + rjc ± yf{Ar 2 K cos 2 fo - A 2 cot 2 <p 3 ). (16) 
19 



Here we obtain tk = 0.201 ± 0.037 from the PQCD analysis 
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and A, 



o 



-0.11 ± 0.065 



by combining recent measurements on CP asymmetry of B° — » K + tt : A cp (B° — » K + tt 



-11.7 ±2.8 ±0.7% 



54, 5j]] with present world averaged value of i?^ = 1.10 ± 0. 15 . 




1.25 



0.95 



FIG. 12: Plot of R K versus with r K = 0.164,0.201 and 0.238. 



As shown in Fig. 12, we can constrain </> 3 with 1 a range of World Averaged Rk as follows: 

• For cos5 >0,r K = 0.164: we can exclude 0° < ^ 3 < 6° and 24° < 4> 3 < 75°. 

• For cos5 >0,r K = 0.201: we can exclude 0° < (f) 3 < 6° and 27° < 3 < 75°. 

• For cos5 > 0, r K = 0.238: we can exclude 0° < 3 < 6° and 34° < 3 < 75°. 

• For cos5 < 0, rx = 0.164: we can exclude 0° < <ps < 6°. 

• For cos5 <0,r K = 0.201: we can exclude 0° < 3 < 6° and 35° < 4> 3 < 51°. 

• For cos5 <0,r K = 0.238: we can exclude 0° < 3 < 6° and 24° < 3 < 62°. 

According to the table 2, since we obtain 5 pi = 157° and <5^/ = 1.4°, the value of cos5 becomes 
negative, —0.91. Therefore the maximum value of the constraint bound for the 4> 3 is strongly 
depend on the uncertainty of \V u b\. When we take the central value of rx = 0.201, c/> 3 is allowed 
within the ranges of 51° < c/> 3 < 129°, which is consistent with the results by the model-independent 
CKM-fit in the (p, rj) plane. 
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VIII. RADIATIVE B-DECAYS (B -» (K*/p/u)i): 



Radiative B-meson decays can provide the most reliable window to understand the framework 
of the Standard Model(SM) and to look for New Physics beyond SM by using the rich sample of 
B-decays. 

In contrast to the inclusive radiative B-decays, exclusive processes such as B — > K*^f are much 



easier to measure in the experiment with a good precision 



6C|. 



Decay Modes 


CLEO 


BaBar 


Belle 


BV{B -> K*°j) (10~ 5 ) 


4.55 ±0.70 ±0.34 


4.23 ±0.40 ±0.22 


4.09 ±0.21 ±0.19 


BV{B -» K* ± 7)(10" 5 ) 


3.76 ±0.86 ±0.28 


3.83 ±0.62 ±0.22 


4.40 ±0.33 ±0.24 


BV{B -> p° 7 ) (HT 6 ) 


< 17 


< 1.2 


< 2.6 


BV(S -► p+7) (10~ 6 ) 


< 13 


< 2.1 


< 2.7 


£V(5 -» o>7) (10~ 6 ) 




< 1.0 


< 4.4 


.4cp(B -» if* 7 ) (%) 


8± 13 ±3 


-3.5 ±9.4 ±2.2 


-6.1 ±5.9 ±1.8 


A C p(B -> if*+ 7 ) (%) 






±5.3 ±8.3 ±1.6 



TABLE VII: Experimental measurements of the averaged branching ratios and CP- violating asym- 
metries of the exclusive B — > V7 decays for V = K*, p and oj. 



(a-l) 



t-K- 



(a-2J 





FIG. 13: Feynman diagrams of the magnetic penguin(a), chromomagnetic penguin(b), annihila- 
tion^) and 02-penguin contributions for B — ► V7 decays 



The main short-distance (SD) contribution to the B — > if* 7 decay rate involves the matrix 
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element 



< K*j\0 7 \B >= |^(-2i)e$ < K^sa^il - j 5 )b\B(p) >, (17) 



which is parameterized in terms of two invariant form fectors as 

<K*(P 3 ,e 3 )\sa^(l- 75 )b\B(P)> = [e 3jM (g ■ P) - P^q ■ e 3 )] ■ 2T 2 (q 2 ) 

+ie^af3e u 3 P a q fS ■ 2Ti(g 2 ). (18) 

Here P and P3 = P — q are the B-meson and K* meson momentum, respectively and £3 is the 
polarization vector of the K* meson. For the real photon emission process the two form factors 
coincide, Ti(0) = 72(0) = T(0). This form factor can be calculable in the kx factorization method 
including the sudakov suppression factor and the threshold resummation effects. As discussed in 
ref 3) we obtain T(0) = 0.28±0.02 for B -> K*j which is far away from the QCD result 0.38±0.06 
by using the light-cone QCD sum rule [olj |. however in accordance with the preliminary result of 
Lattice QCD, 0.25 ± O.OeQ- 

Even though theoretical predictions for the exclusive decays always has large model dependent 
hadronic uncertainties, such uncertainties can be cancelled in the searching of the CP-asymmetry 
and the isospin breaking effect. 

Including all possible contributions from Oj^, 0s g , 02-penguin and annihilation in Figure. 13, we 



obtain the Branching ratios 



63, 



65J: 



• Br(B° -» K°* 7 ) = (3.5±J;£) x 10~ 5 ; Br(B+ -> K+* 7 ) = (3.4±J;|) x 10~ 5 , 

• Br(B° -> p°7) = (0.95 ± 0.14) x 10~ 6 ; Br{B + -> p+j) = (1.63 ± 0.40) x 10~ 6 , 
and the CP- Asymmetry : 

• Acp(B° -> K°* 7 ) = (0.391°;^)% Acp(B+ K+*~/) = (0.62 ± 0.13)% 

The small difference in the branching fraction between _K"°* 7 and K + *^ can be detected as the 
isopsin symmetry breaking which tells us the sign of the combination of the Wilson coefficients, 
Cq/c-j. We obtain 

ri T Br ( B -» K°*~/) - Br (B -» K*~-y) , ,^^ 

where r/ T = t b +/t b o. The first error term comes from the uncertainty of shape parameter of the 
B-meson wave function (0.36 < lub < 0.44) in charm penguin contribution and the second term 
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is origined from the uncertainty of rj T . By using the world averaged value of measurement and 
t b +/t b o = 1.083 ± 0.017, we find numerically that A _(K*7) ex ' p = (3.9 ± 4.8)%. In PQCD large 
isospin symmetry breaking in B — > K*j system cannot be expected. 

IX. SUMMARY AND OUTLOOK 

In this paper I have summarized ingredients of fc^-factorization approach and some important 
theoretical predictions by comparing exparimental data, which is based on my previous works [lI 

y q. 

The PQCD factorization approach provides a useful theoretical framework for a 
systematic analysis on non-leptonic two-body B-meson decays including radiative decays. Our 
results are in a good agreement with experimental data. Specially PQCD predicted large direct 
CP asymmetries in B — > 7r + 7r~, K + it~ decays, which will be a crucial touch stone to distinguish 
our approach from others in future precise measurement. Recently the measurement of the direct 
CP asymmetry in B -» K+tt*, A cp (K + n-) = -12 ± 3% is in accordance with our prediction. 

We discussed the method to determine weak phases (f>2 within the PQCD approach through 
Time-dependent asymmetries in B° — ► 7r + 7r~. We get interesting bounds on 55° < 4>2 < 100° with 
90% C.L. of the recent averaged measurements. 
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